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Simultaneous diffusive and inertial motion of Brownian particles in laminar Couette flow is in- 
vestigated via Lagrangian and Eulerian descriptions to determine the effect of particle inertia on 
diffusive transport in the long-time limit. The classical fluctuation dissipation theorem is used to 
calculate the amplitude of random-force correlations, thereby neglecting corrections of the order of 
the molecular relaxation time to the inverse shear rate. In the diffusive limit (time much greater 
than the particle relaxation time) the fluctuating particle-velocity autocorrelations functions are 
found to be stationary in time, the correlation in the streamwise direction being an exponential 
multiplied by an algebraic function and the cross correlation non-symmetric in the time difference. 
The analytic, non-perturbative, evaluation of the particle-phase total pressure, which is calculated 
to be second order in the Stokes number (a dimensionless measure of particle inertia) , shows that the 
particle phase behaves as a non-Newtonian fluid. The generalized Smoluchowski convective-diffusion 
equation, determined analytically from a combination of the particle-phase pressure tensor and the 
inertial acceleration term, contains a shear-dependent cross derivative term and an additional term 
along the streamwise direction, quadratic in the particle Stokes number. The long-time diffusion 
coefficients associated with the particle flux relative to the carrier flow are found to depend on par- 
ticle inertia such that the streamwise diffusion coefficient becomes negative with increasing Stokes 
number, whereas one of the cross coefficients is always negative. The total diffusion coefficients 
measuring the rate of change of particle mean square displacement are always positive as expected 
from general stability arguments. 

PACS numbers: 05.40.Jc, 47.55.Kf, 82.70.-y, 05.20.Jj 



I. INTRODUCTION 

Combined inertial and diffusive motion of Brownian particles in a flowing fluid is important in a number of aerosol 
processes, including filtration, aerosol sampling, deposition in bends, and particulate deposition in the human respi- 
ratory tract. The limit of negligible inertial effects, where the particles follow closely the motion of the fluid, has been 
extensively studied in sheared colloidal suspensions 0,0- These two-phase systems, which consist of a dispersed par- 
ticulate phase and a continuous fluid phase (gas for aerosols), have long been of interest for their important industrial 
and engineering applications. Pure diffusive particle motion is usually described by a convective gradient-diffusion 
equation from which inertial effects are absent. The limit of inertial transport, where diffusion is neglected and 
particle trajectories deviate significantly from the fluid streamlines, is most conveniently described in terms of the 
particle equations of motion in a Lagrangian formulation. In the transition regime between the diffusion limit and 
the inertia-dominated limit the two particle-transport mechanisms have to be considered simultaneously. 

The effect of particle inertia on the diffusive motion of non-interacting Brownian particles in nonuniform fluids 
has been examined via numerous approaches. Continuum descriptions in terms of mass and momentum conservation 
equations (macroscopic "hydrodynamic" equations) require constitutive relations for the particle-phase total pressure 
tensor. A frequently made approximation is to consider the Brownian particles as an "ideal gas" , thereby using a 
phenomenological expression for the particle-phase pressure tensor 0, Q . 

Mesoscopic descriptions of Brownian motion involve stochastic particle equations of motion and the associated 
Fokker-Planck equation, as, for example, in Subramanian and Brady [f| where a multiple scale analysis of the Fokker- 
■ Planck equation in a simple shear flow is presented. In mesoscopic descriptions in terms of Langevin equations the 
fluctuation dissipation theorem (FDT) is an essential ingredient of the calculation. Santamaria-Holek et al. used 
an alternative approach by considering the motion of a Brownian particle in an external flow field as an example of a 
driven, far-from-equilibrium system. They used mesoscopic nonequilibrium thermodynamics, an approach that does 
not require the specification of the stochastic properties of the random force, to obtain the Fokker-Planck equation 
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for the nonequilibrium Brownian particle distribution function in a simple shear flow. They found that the diffusion 
tensor in the Fokker-Planck equation depends on the shear rate, concluding that fluctuations about the nonequilibrium 
steady state lead to a violation of the classical (equilibrium) FDT. 

Kinetic theory has also been used to obtain the Fokker-Planck equation for the motion of Brownian particles in 
rarefied nonuniform gases. Fernandez de la Mora and Mercer [3 expanded the Boltzmann collision operator in the 
ratio of the light-gas molecular mass to the Brownian-particle mass; they approximated the light-gas distribution 
function by the first two terms in the Chapman-Enskog expansion to derive a Fokker-Planck equation whose diffusion 
tensor was found independent of light-gas velocity gradients p|. However, Rodrigez et al. |8|, using a similar low 
mass-ratio expansion of the collision operator for Maxwellian fluid molecules, obtained a Fokker-Planck equation in a 
uniform shear flow with a diffusion tensor that depended on the irreversible fluid stress tensor. 

Herein, we consider a dilute suspension of non-interacting Brownian particles in a two-dimensional, simple shear 
(laminar, plane Couctte flow). Incrtial Brownian particles arc considered, namely particles whose Stokes number (the 
ratio of the particle relaxation time to the inverse shear rate, or more properly strain rate, a dimensionless number 
that describes particle response to spatial changes in the carrier flow velocity) is finite (at least of order unity). The 
limit of small Stokes number corresponds to pure diffusive motion, whereas for large Stokes numbers inertial transport 
dominates. A mesoscopic approach is adopted, whereby both Lagrangian and Eulerian descriptions are presented to 
investigate the effect of particle inertia on the convective-diffusion equation (equivalently, the Smoluchowski equation) 
in the long-time, diffusive limit. The classical FDT is used since we argue that the modification of the classical FDT 
derived in Refs. [y, |8j is of the order of the ratio of the molecular (fluid) relaxation time (a molecular time scale) to 
the inverse shear rate (a macroscopic time scale), and thus negligible for typical shear rates. 

A linear flow field was chosen as the underlying carrier-gas velocity field because the stochastic particle equations 
of motion and the associated Fokker-Planck equation can be solved exactly for linear flow fields. Thus, the closure 
problem associated with the continuum equations for the dispersed phase is avoided, and analytic, non-perturbative 
expressions may be derived in the long-time limit (for example, for the mean particle velocity). The specific case of a 
simple shear is investigated because Brownian motion in a simple shear has been studied extensively, see, for example, 
Refs. HH 1@. These analyses have been limited in either perturbative calculations 0, @, or in calculations 
of the long-time behavior of equal-time correlation functions without considerations of the associated convective- 
diffusion equation |9(. The long-time, analytic expressions obtained herein extend these previous analyses (under 
well-specified approximations): for example, the non-perturbative, generalized Smoluchowski equation for diffusing, 
inertial Brownian particles in a simple shear is derived, see Ref. Q, and analytic expressions for the particle-phase 
total pressure tensor and the diffusion tensor are obtained, see Ref. 0. Moreover, the methodology presented for a 
simple shear may be used to investigate Brownian motion of finite-inertia particle in any other linear two-dimensional 
flow, be it a symmetric or antisymmetric (rotational) shear. 

Our approach has much in common with the so-called PDF approach used to describe particle dispersion in in- 
homogeneous turbulent flows [10j. Moreover, a Langevin equation formally equivalent to the one used in this work 
has also been used to model turbulence pj| by considering an analogy between the action of the dissipating scales of 
turbulence and that of the molecular white noise driving force in Brownian motion. Hence, the results obtained herein 
apply mutatis mutandis for turbulent dispersion in a simple shear with the proviso that FDT is not a property of 
the turbulent motion. Similarly, a Langevin equation equivalent to the one presented in this work for the fluctuating 
Brownian particle velocities has also been used to describe the effect of shear on fluid velocity fluctuations In 
the identification of similarities and differences between Brownian-particle motion and fluid-point motion the Stokes 
number used in this work corresponds to the dimensionless ratio of the fluid time scale to the inverse shear rate. 

The Lagrangian description is presented in Sec. ITU where the Brownian-particle velocity autocorrelation functions 
are derived. Section lllll contains the Eulerian description and the analytic solution of the Fokker-Planck equation. 
The analytic expression for the total particle-phase pressure tensor as a function of the Stokes number is derived. The 
generalized Smoluchowski equation is calculated in Sec. IIVI as well as the Green-Kubo expressions for the diffusion 
coefficients. The conclusions are summarized in Sec. |3 whereas technical details are presented in the Appendix. 



II. LAGRANGIAN DESCRIPTION 



Consider, thus, a Brownian particle of mass m in a two-dimensional unbounded laminar, plane Couette shear flow, 
namely a simple shear flow with the fluid velocity u along the y-direction, u — ax; the shear rate [2^| tensor a has 
only one non-zero element, a xy = a. The general equation of motion for a small rigid sphere in a nonuniform flow [l3| 
simplifies considerably for most incompressible gas-particle systems because particle density is much greater than fluid 
density. Accordingly, the pressure gradient force, virtual mass, Basset history integral, and Faxen's modification to 
Stokes' drag may be neglected. Moreover, the effect of gravity will be neglected since gravitational settling becomes 
significant only for very large particles (r p » 50 fim). The Saffman lift force will also be neglected, it being negligible 
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with respect to Stokes drag for small-diameter, low-inertia particles. Therefore, the particle equations of motion in a 
Lagrangian description become 

du 

— = [3{ci X -v) + f{t) , (1) 

where the time derivative is a total derivative following the moving Brownian particle with v(t) the particle velocity, 
x(t) the particle position, and f(t) the random force per unit particle mass. As argued, the friction force is assumed 
to be the Stokes drag on the particle. Hence, the friction coefficient j3 is the inverse particle relaxation time, f3 — 
l/r p = 9/jf /(2p p r p ), the particle relaxation time being r p , the particle material density being p p , the particle radius 
r p , and the fluid dynamic viscosity /i/ . As in the Langevin description of Brownian motion in a quiescent fluid the 
random force will be taken to be white in time 



(fi(t)f j (t'))=q5 ij S(t-t') (i,j = x,y) , 



(2) 



with zero mean (fi(t)} = and of an unspecified, at the moment, strength q. Angular brackets (•} denote an ensemble 
average over all particle trajectories. 

The Langevin equations are linear: hence, the particle velocity and position may be solved formally as functionals 
of the random force. The formal solution for Brownian particles injected at the origin with zero initial velocity (the 
choice of the initial conditions is not important since we are interested in the long-time behavior) are 



v x (t) = < 

Vy{t) 



dhe^ [(3ay(h) + f x (h)} 



x(t) = I dtx v x (tx) , 



dtl Vy{t\ 



(3a) 
(3b) 



The formal solution allows the analytic evaluation of ensemble averages of products of particle position and velocity, 
two-point correlation functions, in terms of the random-force strength q. Their evaluation requires the calculation of 
time integrals whose integrands include the causal correlation 



(y(t)fy(t')) = 0(t-f)| 



1 - e 



-0(t-t') 



(4) 



The Heaviside 9 function in the previous equation arises naturally via the explicit evaluation of the time-dependent 
correlation function using Eas. 13 bl it ensures causality. 

In the diffusive limit (t ^> /3 _1 ) exponential terms in the equal-time, two-point correlation functions may be 
neglected. If only polynomial terms in time are kept the correlations evaluate to 
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(v(t)v(t)) = 
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(5c) 



In a quiescent fluid the fluctuation strength is specified by invoking the fluctuation-dissipation theorem. Its classical 
form may be expressed as (see, for example, Reeks [l4|) 

/>OC 
777 / 

The use of the classical FDT in sheared systems has been questioned. Whereas its classical form as shown in Eq. © 
has been used in the past 0,00]! it has also been argued 0] that the fluctuation strength may be determined by the 
ad hoc requirement that energy equipartition hold in the local, co-moving reference frame. Energy equipartition in the 
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local reference frame implies local equilibrium and hence the Brownian velocity distribution function becomes locally 
Maxwellian. As we show in Sec. 11111 if the Brownian velocity distribution function is a local Maxwellian distribution 
then the Brownian particles behave as an ideal "particle" gas with no shear stresses, nor shear viscosity. 

More recently, Santamaria-Holek et al. derived a Fokker-Planck equation for the nonequilibrium distribution 
function of Brownian particles in stationary flow. They avoided the use of a Langevin equation (and the associated 
problem of the specification of the statistics of the random force) by using arguments based on mesoscopic nonequi- 
librium thermodynamics. They found that the flow modifies the diffusion tensor in the Fokker-Planck equation by 
a term proportional to the imposed velocity gradient. They concluded that due to this additional term the Fluc- 
tuation Dissipation Theorem in its classical form does not hold for fluctuations about the nonequilibrium steady 
state. However, they did not provide an estimate of the relative magnitude of the nonequilibrium correction to the 
classical FDT |29j . Such an estimate may be deduced from earlier kinetic theory calculations of particle motion in a 
nonuniform light gas. Rodriguez et al. Q used Fokker-Planck and Langevin descriptions of fluctuations in uniform 
shear flow to conclude that the diffusion tensor in the Fokker-Planck equation, and the corresponding properties 
of the random force in the Langevin description, are modified by a term proportional to the traceless part of the 
fluid pressure tensor. Thus, the nonequilibrium modification of the FDT implies that in a Langevin description the 
stochastic properties of the random force are modified by the shear flow. However, in a Langevin description the 
time scale of the Brownian white noise driving force is considered much shorter than the time scale of the imposed 
flow suggesting that the nonequilibrium correction would be of the order of the ratio of the molecular relaxation time 
to the time scale of the imposed shear 30]. In fact, the nonequilibrium correction derived by Rodriguez et al. @ 
can be shown to be proportional to the ratio of the fluid molecular relaxation time to the time scale of the imposed 
shear (a ratio of a microscopic to a macroscopic time scale) by expressing the (dimensionless) correction in terms of 
the velocity gradient (shear rate) and the ratio of the fluid viscosity to the fluid pressure. Similarly, Fernandez de la 
Mora and Mercer Q used a Chapman-Enskog expansion of the light-gas velocity distribution function to show that 
the nonequilibrium modification of the Fokker-Planck equation is of the order of the ratio of the light-gas relaxation 
time (rf ~ fJ-f/pf) to the macroscopic fluid deceleration time (the inverse shear rate). Thus, it becomes of the order 
of the light-gas Knudsen number and therefore it may be neglected, only becoming important for rarefied gases 0. 
Since in this work we are interested in cases where time scales are clearly separated, in what follows we shall use the 
classical form of the FDT that does not introduce additional random-force correlations. 

The classical form of the FDT and Eq. specify the fluctuation amplitude to be 

q_ = k B T 

2/3 m U 

With this identification some of Eqs. (J5J have been previously reported |{| , whereas the full time-dependent correlations 
for delta- function initial conditions are given in Ref. |£g. Particle velocity correlations along the y direction (in the 
long-time limit), being independent of the shear flow, satisfy the normal equipartition theorem. This is expected 
since the shear flow is only along the ^-direction, and the (x, y) components of the random force have been assumed 
uncorrelated. Moreover, a kinetic temperature, as opposed to the thermodynamic temperature T in Eqs. i|5c|) . may be 
defined by relating the average particle kinetic energy to fcsTkin- Then, the kinetic temperature along the direction of 
the shear becomes time-dependent and a function of particle mass (through the dependence on /3) and of properties 
of the fluid flow (the local fluid gradient a). 

The long-time limit of the equal time correlation (y 2 (t)) may be used to obtain the Stokes-Einstein expression for 
the diffusion coefficient, 

D °=2jyw = ik ■ (8) 

Note that the diffusion coefficient is first order in /3 _1 , a result that will be used later in Sec. I1VI 

The shear flow also modifies the time-dependent particle-velocity autocorrelation functions. Ensemble averages 
of the formal solutions of the equations of motions Eqs. ((3J) give in the diffusive limit (t 3> /? _1 ), neglecting terms 
involving powers of exp(— j3t), 

, , fifl^ \e- f3T + ± St 2 [4/3i - 8 - e-^ {fir + 3)11 for r > , 
(vJt + r)v x (t)) = < , m T L 2 Y > * , sn ~ (9a) 

V V ' y " \ e^ T + ±St 2 4/3 (t + r -8 + e^ T /3r - 3 for r<0 , V ' 
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(Vy(t + T)Vy(t)} 



-ft\r\ Vt 



(9c) 



Hence, the combined effects of particie inertia and shear flow modify both the amplitude of the autocorrelation 
functions and their time dependence. The dependence of the autocorrelation functions on particle inertia and shear 
rate has been expressed in terms of the Stokes number, St = a/ f3. The shear flow not only breaks spatial symmetry but 
also (macroscopic) time reversibility and stationarity: the particle- velocity autocorrelation function in the streamwise 
direction is nonstationary. Of course, the velocity correlations perpendicular to the shear decay exponentially in 
time as expected for a Gaussian, stationary Markov process, an Ornstein-Uhlenbeck process. The shear-induced 
modifications of the autocorrelation functions become more transparent if the contribution of the underlying shear 
flow is subtracted. Specifically, the autocorrelation functions of the fluctuating particle velocities with respect to the 
shear flow, v'^(t) = v x (t) — ay(t), become (in the long-time limit) 



«(rK'(0)> = 
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(10b) 



Thus, the autocorrelation functions of the fluctuating particle velocities are stationary, the velocity correlation along 
the shear is symmetric in the time difference r, but the cross correlation is non-symmetric in r. Since in the local 
reference frame time stationarity is recovered the dependence of the autocorrelation functions on t has been dropped: 
note, however, that these expressions are valid only in the long-time limit (f3t ^> 1). The time decay of the velocity 
correlation function along the flow direction is not a pure exponential, but it is modified by an algebraic prefactor; 
hence, the underlying stochastic process v"(t) is not an Ornstein-Uhlenbeck process. Furthermore, since the cross 
correlation is stationary (v"(t)v v (0)) = (v"(0)v v (—t)), but since it is not symmetric (v"(t)v v (Q)) ^ (v^(0)v v (t)). 
The time asymmetry of the cross correlation implies that for a negative time difference r < the correlation decays 
exponentially, whereas for r > the correlation initially increases at short relative times to decrease exponentially at 
longer time. The maximum occurs at j3r — 1/2. 

These qualitative observations are summarized in Fig. ^ where the three correlation functions are compared. The 
normalized (at the origin) autocorrelation functions are plotted, Cy (r) = (v" {t)v" {fy) / {v" {fyv'j {$)) as functions of 
time rendered dimcnsionless by a (3 scaling. The autocorrelation along the streamwise direction C xx (t) has been 
parametrized by three values of Stokes number; the other two normalized correlations, C xy (r) and C yy (r) are in- 
dependent of the Stokes number. Note that with increasing Stokes number the cusp at the origin, characteristic of 
the short-time behavior of C TO (r), becomes rounded. The St = 10 curve is also the value of C xx (t) in the limit the 
St — > oo. The velocity autocorrelation functions will be reconsidered in Sec. IIVI since they are related to particle 
diffusion coefficients through Green-Kubo relations. 

Reeks m an analysis of particle turbulent dispersion in a shear flow reports equations similar to Eqs. (|10fl . It can 
be shown that the results of Ref. 01 f° r the case of fluid-point dispersion lead to velocity autocorrelation functions 
identical to those presented above. Similarly, Eckhardt and Pandit 12] obtained results formally equivalent to ours 
for the effect of shear on fluid velocity fluctuations. Note, however, that for the system under consideration the effect 
of the shear flow and particle inertia are described in terms of the particle Stokes number: hence, for a given shear 
rate changes of the Stokes number imply changes of the particle relaxation time (i.e. of the particle's inertia). 



III. EULERIAN DESCRIPTION 



The Eulerian description of the motion of N, independent and identical Brownian particles in terms of the aver- 
age phase space density P{x, v; t) is an alternative to the Lagrangian description in terms of stochastic differential 
equations. The Fokker-Planck equation associated with the particle equations of motion Eqs. Q may be derived in 
numerous ways, including the elegant functional derivation of Refs. |!8lll9| . It becomes 



m dv 2 



(11) 



where the fluctuation amplitude, defined in Eq. (JJJ, was used. The probability density function P is normalized on 
the total number of particles. As mentioned before, Eq. (|llf> differs from the Fokker-Planck equation derived by 
Santamaria-Holek et al. B in the diffusive term since we used the classical FDT. 
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FIG. 1: Time dependence of the fluctuating particle-velocity autocorrelation functions in the diffusive limit, parametrized by 
the Stokes number: solid lines are C xx (t), the dashed line C xv (t), and the dash-dot line is C vv {t). 



The probability density function P is used to define average quantities. In particular, the mean, local particle 
velocity, a quantity defined as in kinetic theory of gases |2l| is 

p(x;t) ~v(x;t) = m J dv P(x,v;t)v = m{nv) , (12) 

with n the "instantaneous" particle number density and p(x;t) the local, average particle mass density, 

p(x;t) = m J dvP(x,v;t) . (13) 

The fluctuating component of the particle velocity field, defined with respect to the mean particle velocity [3lj , will 
be denoted by v' = v — v. 

These local averages, denoted by an overbar and known as density-weighted averages in two-fluid descriptions of 
dispersed flows, are averages with respect to a local, properly normalized velocity probability density distribution 

P(x,v;t) 



ip(x,v;t) 



J dv P(x, v; t) 



(14) 



The integral of if) over particle velocities is, thus, unity. 

The appropriate equations to describe particle dispersion, also known as "continuum" or "hydrodynamic" equations, 
are obtained by taking velocity moments of the Fokker-Planck equation. In particular, the particle mass conservation 
equation is obtained by multiplying Eq. (|11J) by m, integrating over particle velocities and using Eq. (H3J. Multiplying 
the Fokker-Planck equation by mv' and realizing that m (nv'v') = pv'v' [see the definition Eq. I)14[l] gives the 
momentum conservation equation. Accordingly, they become (see also |l8|) 

dp 
~dl 



-pV ■ v , 



(15a) 
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where the total derivative d/dt — d/dt + v ■ V describes changes with respect to the mean particle velocity. Equa- 
tions l|15f) are a special case for a white noise random force of the general continuum equations for the dispersed phase 
in inhomogeneous flows derived by Reeks [T(i| . 

The momentum equation Eq. (|1 5 b>|) may be compared to the momentum equation that expresses momentum balance 
in the particulate phase [22, 

^.=(3(u-v)--V-P p , (16) 
dt p 



where P p is the unspecified, particle-phase total pressure tensor. Comparison yields P p = pv'v' . It is easy to show 
[via Eq. (|14p] that the expression for the particle-phase pressure tensor is identical to the pressure tensor as defined 
in kinetic theory |2l| or in mesoscopic nonequilibrium thermodynamics 



J dv(v - v)(v -v)P(x,v;t) . (17) 



It is apparent from the previous discussion that the continuum equations and the identification of the particle-phase 
pressure tensor with the particle covariances only require that the random force be white in time. Thus, the previous 
expressions are valid for a general flow field and not only for a linear flow field (as long as the random force is white). 
In what follows we restrict the calculation to a linear flow field, eventually evaluating average quantities (for example, 
the mean particle velocity and the particle-phase pressure tensor) for the specific case of a simple shear flow, as 
discussed and justified in the Introduction. 

For a linear flow field the Fokker-Planck equation defines a Gaussian process for [x,v] and thus it has an analytic 
solution, see, for example, Refs. [l8l , l2Cj . The Gaussian solution may be used to evaluate explicitly the density-weighted 
ensemble averages; Swailes and Dlrbjshire M report the analytic expressions. For Brownian particles injected at the 
origin of the co-ordinate system with zero initial velocity the spatially-dependent particle concentration is (Ref. [2£j 
presents the analytic solutions for arbitrary initial conditions) 



/ \ 1 
p(x;t) = 7T7 ex P 

2tt [det((a%c})] 1/2 



1 T / \—l 

—x ■ (xx) ■ X 

2 v 1 



(18) 



the mean particle velocity 

v(x;t) — (vx) ■ {xx)^ 1 ■ x , (19) 

and the particle-velocity covariances 



v'v'(t) = (vv) - (vx) ■ (xx)- 1 ■ (xv) , (20) 

where the superscript in Eq. (|18fl denotes transpose. Thus, for the Gaussian process the mean particle velocity is 
linear in x, and the particle covariances are spatially independent. Analytic expressions for the long-time behavior of 
the mean-particle velocity and its approach to the steady-state value (the carrier flow velocity) are presented in the 
Appendix. For completeness note that explicit evaluation of the necessary correlation functions in the long-time limit 
shows that lim^oo v'v' = ((v — u) 2 ) since lim t _ i00 (v) = u. 

The long-time particle-phase pressure tensor may then be evaluated using Eqs. JSJ and l|20() : expressed as a function 
of the Stokes number, it becomes 



P P = P 



1 



St 2 \ « St /o l\ St 2 (\ 



2 V 1 0/ 4 \0 -1 



(21) 



The particle pressure tensor has been decomposed into a form reminiscent of the pressure tensor of a simple liquid 
undergoing laminar plane Couette flow 0- The first term is isotropic and proportional to the identity tensor; the 
other two terms constitute what has been called [24[ the "friction pressure tensor" [33] • The first part of the isotropic 
term gives the ideal gas pressure of a collection of Brownian particles ("Brownian particle gas"), a result that has 
been previously postulated on phenomenological arguments [3|, |j] ■ The ideal gas pressure is modified by a correction 
dependent on the Stokes number. The other terms, usually neglected, arise from the particle viscous stresses. The 
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second term is proportional to the symmetric rate-of-strain fluid tensor, 2 *e = VS + fv^j • Thus, the (negative) 
proportionality constant defines the shear viscosity of the particle phase to be 

r,t« ai = \D p . (22) 

Hence, the Brownian shear viscosity, being independent of the shear rate, is a conventional Newtonian viscosity. An 
alternative expression for Eq. (|22|l is that the particle Schmidt number, the ratio of momentum diffusivity to mass 
diffusivity, is Sc p — (rjp ax /p)/Do = 1/2, as is the case for particle motion in a turbulent flow in which the particle 
response time is much larger than the turbulent time scale 17] . The third term, which is absent in Newtonian fluids, 
shows that the particle phase exhibits non-Newtonian rheology, as has been remarked for sheared simple liquids |2U24|. 
Note that the last term, as well as the correction to the ideal gas pressure, is second order in the Stokes number. Hence 
to leading order in St the viscous part of the particle pressure tensor is traceless, and the particle phase behaves as 
a Newtonian fluid. Non-Newtonian behavior becomes evident only to second order. Moreover, according to Eq. 121|) 

[and Eq. ||SJ|] in the absence of diffusion, for sufficiently massive particles, P p — 0. 

As remarked in Sec. [H]if local equilibrium is assumed, as for example in Ref. |l6j . ((v — u)(v — it)) = fc^Tl/m, 
the Brownian particle distribution becomes locally Maxwcllian. Then, from Eqs. (|5c|l the fluctuation strengths 
can be calculated to form a symmetric tensor dependent on the shear rate. As extensively argued in Sec. [H] the 
dependence of the fluctuation strengths on the shear rate implies that the externally imposed shear modifies the 
properties of the random force, a condition that is unlikely to hold when the relevant time scales (inverse shear rate 
and molecular relaxation time scale) are clearly separated. Moreover, for a local equilibrium assumption a similar 
calculation shows that the particle-phase pressure tensor (in the long-time limit) contracts to the ideal gas result, 

P p eq ' = p(x;t)kBTl/m, i.e. there is no viscous shearing. 

Equation (|21|l along with Eqs. (|15|l are our main result. They are analytic, non-perturbative expressions that 
describe the long-time behavior of the dispersed phase: they incorporate the effects of the flow field, and the combined 
effects of particle diffusion and particle inertia on the particle-phase total pressure tensor. However, the momentum 
equation is coupled to the mass conservation equation. The two (steady-state) equations are, usually, decoupled by 
performing a low Stokes-number expansion of the momentum equation to obtain a perturbative expression for the 
mean particle velocity field in terms of the fluid velocity and its gradients 0]. Substitution of the resulting mean 
particle velocity into the particle mass conservation equation then gives the associated convective-diffusion equation [j| . 
For the case of a simple shear this procedure leads to a generalized Smoluchowski equation [since the local average 
density is defined as an integral of the average phase space density P, Eq. 1)13(1]. as we show in the following Section. 



IV. SMOLUCHOWSKI EQUATION 

The momentum equation (|15b|) . an equation valid for a general flow field and a white random force, may be 
re-arranged to obtain an explicit expression for the total time derivative of the moan particle-velocity field (also 
referred to as the inertial acceleration term). For the specific case of a simple shear differentiation of the analytic, 
time-dependent expression for the particle concentration Eq. l(T%|) leads to 

px = -(xx}-^ , (23) 

which substituted into the equation for the mean particle velocity Eq. 1(19(1 expresses the particle mean velocity in 
terms of the density gradient 

v = -(vx)--^- . (24) 
p ox 

Moreover, a variant of Eq. 1(23(1 . obtained by multiplying it by the shear-rate tensor S, determines the carrier flow 
velocity 

pu = -{ux}-^ . (25) 

Equations ((24(1 and 1(25(1 . derived from the time-dependent solution for the local particle density, arc generally valid 
for a simple shear and not only in the long-time limit. Their substitution into the momentum equation, Eq. I(15b(l . 
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along with the realization that for a linear flow field (and, in particular, a simple shear) the velocity covariances are 
spatially independent, leads to 

4_[«(.-.>. > -roi.g*fl.& . (26, 

Hence, the inertial acceleration term in a simple shear may be explicitly evaluated: specifically, Eqs. JHJ along 

with Eq. J7J and Eq. (|2U[1 provide analytic expressions for the necessary correlations (in the diffusive limit). Their 
substitution into Eq. (|26|l yields 

M=MW-f 'A . (27) 

Thus, the proportionality matrix M is of the order of the Stokes number; note, however, that for a simple shear the 
convective derivatives of the flow field [(«• V)it] vanish. In this respect the simple shear differs from the symmetric and 
antisymmetric shear, the other possible representations of a general two-dimensional linear flow field. Equation 1271) 
for the inertial acceleration term and Eq. (|21|) for the particle-phase pressure tensor were derived for the specific 
case of a simple shear. Hence, they may be substituted into the general, particulate-phase, momentum-conservation 
equation Eq. JTBJ, or equivalently into Eq. i|15b|) . rewritten as 

*="-s v -Mf ' (28) 

to obtain the analytic, non-perturbative expression for the mean particle velocity field in the diffusive limit (t 3> /3 _1 ) 

v = u - ((v - u)x) ■ (29) 



= U- D 



p dx 



l 7T ■ < 3 °) 

p OX 



Note that Eq. (|29(l may also be derived by a judicious combination of Eqs. i|24[l and i|25|l as can be seen by expanding 
the ensemble average ((v — u)x). The alternative derivation presented earlier, however, is more instructive because it 
allows the identification of the contribution of the particle-phase pressure tensor and the inertial acceleration to the 
diffusion coefficient, as discussed in what follows. 

To leading order in the particle relaxation time (/3 _1 ) the mean particle velocity is the fluid velocity modified by 
a diffusion term (Fick diffusion) . Fick diffusion stems from the ideal part of the particle-phase pressure tensor. The 
correction to Fick's law arises from the inertial acceleration and the particle viscous stresses. In fact, the contribution 
of the inertial acceleration term in Eq. (|28|l is of the same order as the contribution from the particle-phase viscous 
stresses. Either contribution to the mean-particle velocity is at least second order in /3 _1 [see, also, Eq. (JSJ)]. Hence, 
in a simple shear it is inconsistent to retain particle viscous stresses and neglect inertial acceleration, and vice versa. 

The convective-diffusion equation (Smoluchowski equation) associated with the two coupled Eqs. (|15[) is obtained 
by substituting the mean particle velocity Eq. (|29|1 into the mass conservation equation to obtain 

d ± + V.{pu) = {{v-u)x).^ . (31) 

Equation (|31|l defines the diffusion tensor D = (v"x) that depends on both particle inertia and the shear rate through 
the Stokes number. In the long-time limit it is given by the second term of the RHS of Eq. (|30[1 . Note that the diffusion 
tensor, as well as the matrix expression for the inertial acceleration term Eq. (|27|) . is not symmetric, reflecting the 
symmetry-breaking effect of the imposed shear. Therefore, the long-time generalized Smoluchowski equation becomes 



Since the inertial acceleration term is anti-symmetric it only modifies the diffusion coefficient in the streamwisc 
direction; the modification is second order in the Stokes number. The first order correction through the appearance 
of the cross derivative term arises solely from the viscous part of the particle pressure tensor. The same first order 
correction to the Smoluchowski equation was derived by Subramanian and Brady pj, who, however, did not calculate 
higher order corrections. 
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cT 




St=o/p 

FIG. 2: Dependence of the long-time diffusion coefficients on particle inertia expressed in terms of the Stokes number. The 
components D yx (dash-dot line) and D yy (dotted) are always positive, D xy (dashed) always negative, whereas the sign of D X x 
(solid) changes as a function of the Stokes number. 



The sign of the long-time diffusion coefficient along the flow direction, D xx in Eq. \YVll , depends on the value of the 
shear rate. For large shear rates it becomes negative, the critical Stokes number being St = a/3/2. The long-time 
limit of the off-diagonal diffusion coefficients is such that D yx is always positive, whereas D xy is always negative [See 
Eq. 130fl ]. The other diffusion coefficient, D yy , is always positive. The dependence of the diffusion coefficient on 
particle inertia (via the Stokes number) is shown in Fig. 

The existence of negative diffusion coefficients does not violate any stability conditions since the coefficients that 
appear in the convective-diffusion equation refer to the particle diffusive flux with respect to the underlying carrier 
flow. As expected from general stability arguments the total diffusive flux defined with respect to a fixed reference 
frame, and not with respect to the carrier flow, is always positive. The mean particle velocity expressed in terms of 
the density gradient, Eq. 124(1 . defines the total diffusion coefficients, that in the diffusive, long-time limit become 



Dtot = (vx) = D 



1 



2/2 



-aH 

— 





2at 



(33) 



Therefore, the total diffusion tensor is time dependent, non-symmetric, the corresponding matrix positive definite, 
and its components are always positive, thereby assuring the stability of the system. The shear-induced modification 
of the total diffusion coefficients, in particular the quadratic time dependence of the strcamwisc diffusion coefficient, 

has been noted before, for example in Refs. @>|UE|- Note that the symmetric part of D to t gives the total diffusion 
coefficients as determined from the rate of change of the particle mean square displacement. 

The long-time diffusion coefficients may also be obtained from the velocity autocorrelation functions presented in 
Sec. [H] The following Green-Kubo relations, also known as Taylor's diffusion formulae for diffusion via continuous 
movements in theories of turbulent particle dispersion, express the diffusion coefficients in terms of the velocity 
autocorrelation functions 



5(t) = (v"(t)x(t)) = f dh (v' 
Jo 



(*)«(*!)> 



(34) 
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In the diffusive limit they become 



/ 

Jo 



oo 




dti lim [(v"(t)v(t - h))] 



(35) 



expressions that reproduce the diffusion tensor Eq. (|3(JI) [see, also, Eq. (| A. 5|> ] . 

Equation (|35l) . in addition to providing the long-time diffusion coefficients in terms of velocity autocorrelations, 
allows a physical interpretation of the negative diffusion coefficients. Specifically, the long-time limit of D xy is always 
negative because a particle crossing the mean shear flow with positive v y will have negative streamwise fluctuations (as 
the mean flow increases with y). Along the streamwise direction D xx becomes negative for large shear rates because 
streamwise fluctuations v x = v x — ay become negative for large values of the shear rate a (for v x > 0). 

The results presented in this Section may be critically compared to earlier analyses of diffusive motion in nonuniform 
flows. We showed, via the derivation of analytic, long-time expressions that in a simple shear the inertial acceleration 
term (total derivative of the mean particle velocity) and the particle-phase viscous stresses have to be treated consis- 
tently. In particular, we argued that in perturbative evaluations of, for example, the Brownian particle diffusive flux 
in a simple shear, it is inconsistent to retain the particle-phase viscous stresses and neglect inertial acceleration, and 
vice versa. Previous studies differ from ours in the way these two terms (P p and dv/dt) are treated (and to a lesser 
degree in the choice of the underlying flow field). 

Specifically, Fernandez de la Mora and Rosner £ 4] considered the effect of inertia on diffusional deposition for a 
general flow field. They solved the steady-state momentum equation, Eq. (|28(l . to leading order in keeping only 
the ideal-gas part of P p and the leading order contribution from the inertial acceleration term (it • V)it]. For 
a general flow field this is a consistent approximation: higher order corrections would require an expression for the 
particle viscous stresses. For a simple shear their result reduces to ours to the same order [v = u — D \7 log p+0(f3~ 2 )] 
since, as noted earlier, the convective derivatives of the flow field vanish in a simple shear. 

Ramshaw Q also analyzed Brownian motion in a general flow field, but his analysis differs from ours in that 
particles of negligible inertia were considered. He developed a leading order in expansion neglecting both viscous 
stresses and inertial acceleration, but keeping other forces, e.g. thermophoresis. Ramshaw jl^ extended the original 
derivation to include phenomenologically viscous stresses in the particle and mixture (particle and fluid) momentum 
equations retaining however the assumption of negligible particle inertia. He found that inclusion of viscous stresses in 
the mixture equation, and evaluation of the suspension viscosity by Einstein's formula, leads to an additional term in 
the diffusion tensor proportional to the particle-phase volume fraction. Our calculation considers an infinitely dilute 
suspension, and hence this additional term is absent from Eq. (|3U[) . Moreover, his calculation shows that inclusion of 
viscous stresses only in the particle-phase momentum equation does not modify the diffusion coefficient if the particle 
shear viscosity is evaluated as if the particles were an ideal gas (and, hence, the shear viscosity is independent of the 
particle mass density). 

The analysis of Santamaria-Holek et al. , who considered diffusion in a simple shear as in this work, differs from 
ours primarily in the use of the nonequilibrium FDT (as summarized in Sec.[nJ. Moreover, the diffusion tensor was 
evaluated perturbatively in the particle relaxation time (/3 _1 ) retaining the viscous part of the particle-phase pressure 
tensor and neglecting the contribution from the inertial acceleration term. As shown earlier, when the classical form 
of the FDT is used these two terms contribute to the same order (/3 -2 ) in the perturbation expansion. However, 
the resulting, leading order correction to the Smoluchowski equation is identical to ours (up to the nonequilibrium 
modification) because the inertial acceleration (being anti-symmetric) contributes only to the highest order. A similar 
remark applies to the leading order correction to the Smoluchowski equation derived by Subramanian and Brady Q 
(who neglected the nonequilibrium modification, as in this work). 



The primary question addressed in this work is how particle inertia modifies diffusional transport of particles in a 
nonuniform flow. The usual convective-diffusion equation describes diffusional transport in the limit where the particles 
follow the fluid streamlines. It is an Eulerian continuum equation that provides a computationally efficient method to 
calculate particle transport (and deposition); however inertial effects are neglected. Inertial transport, where particle 
trajectories deviate considerably from the fluid streamlines, is best calculated via a Lagrangian description in terms 
of the particle equations of motion. Diffusive particle motion may be incorporated in the Lagrangian formulation 
through the addition of a random force, but the numerical solution of the resulting stochastic differential equations is 
computationally intensive. It is, thus, desirable to derive a continuum equation valid in the transition regime between 
the diffusion limit and the inertia-dominated limit, incorporating both particle-transport mechanisms. 



V. 



CONCLUSIONS 
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We considered the coupled diffusive and inertial motion of non-interacting Brownian particles in a simple inhomo- 
geneous fluid flow, a simple shear flow. The long-time, diffusive, behavior of the system was investigated neglecting 
the short-time regime. Even though the choice of a simple shear as the underlying carrier flow is restrictive [for 
example, the convective derivatives of the flow vanish, (u ■ V)it = 0, a condition that does not hold for a symmetric 
or antisymmetric shear] the choice was motivated by numerous previous investigations of Brownian motion in such a 
flow. More importantly, the Fokker-Planck equation associated with the particle equations of motion in a linear flow 
field is of the linear type 8] , and hence analytically solvable. The solution of the Fokker-Planck equation was used to 
obtain analytic expressions for average particle properties, for example the mean particle velocity and particle- velocity 
correlations. These analytic solutions depend on equal-time ensemble averages that were determined from the formal 
solutions of the Langevin equations for particle motion. The classical form of the Fluctuation Dissipation Theorem 
was used to determine the strength of the random force correlations, neglecting corrections of the order of the ratio 
of the fluid molecular relaxation time to the time scale of the imposed shear. We showed that the long-time, time- 
dependent particle-velocity autocorrelation along the streamwise direction is nonstationary. The fluctuating velocity 
autocorrelations were determined to be stationary in time, but the cross correlation was non-symmetric in the time 
difference, reflecting the combined effect of particle inertia and shear on particle-velocity fluctuations. 

We used the analytic solution of the Fokker-Planck equation, in conjunction with its first two velocity moment 
equations, to obtain [in the diffusive limit (t S> /3 -1 )] the generalized convective-diffusion equation (generalized 
Smoluchowski equation) that incorporates inertial effects on diffusional transport for dilute suspensions. The coupling 
of particle inertia to the fluid flow introduces a shear-dependent, linear in the (particle) Stokes number, cross derivative 
term, and an additional term along the streamwise direction, quadratic in the Stokes number. The associated diffusion 
tensor, thus, depends on the shear rate and particle inertia; the diffusion coefficient along the streamwise direction 
flow was found to become negative for large particle relaxation times (or, equivalently, for large shear rates), whereas 
one of the cross diffusion coefficient was determined to be always negative, the other two being always positive. We 
argued that stability conditions are not violated since the total diffusion coefficients (not those with respect to the 
carrier flow) that measure the rate of change of particle mean square displacement were determined to be always 
positive. 

We showed that in a simple shear, and in the long-time limit, the contribution of the inertial acceleration term to 
the diffusion tensor in the generalized Smoluchowski equation is of the same order in the Stokes number as that of the 
viscous part of the particle-phase pressure tensor. Thus, in perturbative evaluations of, for example, the Brownian 
particle diffusive flux, it is inconsistent to retain the particle-phase viscous stresses and neglect the inertial acceleration 
term, and vice versa. 

As part of the derivation of the Smoluchowski equation we calculated the particle-phase total pressure tensor that 
was determined to be second order in the Stokes number. Similar to the pressure tensor of simple sheared liquids, 
the pressure tensor was decomposed into three parts: a part proportional to the identity tensor that gives the ideal 
pressure of a gas of Brownian particle with an additional term due to particle inertia, and two terms that arise from 
the particle viscous stresses. We found that the particle phase has a conventional shear viscosity, but it behaves as a 
non-Newtonian fluid if second order effects in the Stokes number are considered. 

These results will be extended to other two-dimensional linear flows in future work. 
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APPENDIX: APPROACH TO THE STEADY STATE 



The time-dependent, equal-time correlation functions in conjunction with Eq. I|19|l [or Eq. I|24[l ] may be used to 
calculate the approach of the mean particle velocity to its steady-state value. The required calculations were performed 
with Mathematica [26(. Accordingly, the mean particle velocity approaches the carrier fluid velocity as follows 
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Note that the long-time limit and the vanishing shear-rate limit do not commute since 
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a + 0(a 2 ) 



/3i 4(/3f) 5 



a + 0(a z ) 



(A.2a) 
(A.2b) 



The zeroth order terms of Eqs. I|A.2|) give the mean particle velocity of Brownian particles in a quiescent fluid. 

The time-dependent approach of the mean particle velocity to its steady-state value may be used to estimate the 
relative importance of the two terms in the total derivative of the mean-particle velocity, namely d/dt and (v ■ V). 
Appropriate differentiation of Eqs. (|A.1|I . along with 



<91ogp 1 m 3y 
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St k B T t 2 
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(A.3a) 



^l = -^-+0(t- 2 ) 
dy D t +U[t > 



(A.3b) 



shows that the explicit time derivative vanishes in the long-time limit (to leading order in t 1 ) to give p(v ■ V)t> = 

M ■ dp/dx. Hence, the leading order contribution to the total derivative does not arise from the explicit time 
derivative. 

For completeness we also present the approach to steady-state values of the particle-velocity covarianccs 
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The evaluation of the streamwise diffusion coefficient D xx via the Green-Kubo relations Eq. 135(1 requires the 
nonstationary correlation function 
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